LIMIT THEOREMS FOR SELF-SIMILAR TILINGS 
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Abstract. We study deviation of ergodic averages for dynamical systems 
given by self-similar tilings on the plane and in higher dimensions. The main 
object of our paper is a special family of finitely-additive measures for our 
systems. An asymptotic formula is given for ergodic integrals in terms of these 
finitely-additive measures, and, as a corollary, limit theorems are obtained for 
dynamical systems given by self-similar tilings. 

1. Introduction 

We study the deviation of ergodic averages for certain tiling dynamical systems, 
namely, translation Reactions associated with self-similar tilings. For d = 1 
asymptotic formulas and limit laws for such deviations were obtained in |10| . 
The main novelty of the d > 2 case is the appearance of "boundary effects" which 
result in some new phenomena. 

We assume that the tilings have translationally finite local complexity, are ape- 
riodic, and repetitive. Self-similarity means that there is an expanding similarity 
map (f> : R rf — > such that every "inflated" tile can be subdivided into tiles 
of the original tiling, basically, a Markov property. Given a self-similar tiling, 
we consider its orbit under translations and its closure in the natural "local" 
topology. This is a compact metric space, on which M d acts by translations. The 
resulting dynamical system is known to be minimal and uniquely ergodic. See 
the next section for precise definitions and statements. 

Let S be the substitution matrix, which is primitive, and let 6i,...,9 m be 
its eigenvalues, ordered by their absolute values: 0\ > \02\ > ■ ■ ■ > \9 m \. In 
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order to describe our results, we make a simplifying assumption that there are no 
Jordan blocks associated with eigenvalues of absolute value \02\. The most basic 
class of functions for which we consider the deviation of ergodic averages is the 
collection of characteristic functions for "cylinder sets" of tiles. Averaging can 
be done over balls or cubes of diameter R, or more general increasing families 
of Lipschitz domains. A question arises: how can we estimate from above the 
deviation of the ergodic average from the mean? The answer depends on the 

d-1 d-1 

relation between \02\ and l d , see Corollary |4.51 If 1 6*2 1 < ®\ d > then the deviation 
term is bounded above by Ci?^ -1 , which means that the main contribution comes 

d-1 

from the boundary of the domain. On the other hand, if 1 6*2 1 > ®\ d , then the 
deviation term is bounded above by CR a , where a = ^|2I^J g (^ _ \^d) (if 

d-1 

I $2 1 = Q\ d > then there is a logarithmic correction). These deviation bounds are 
sharp, at least, in the general case. There are related recent results by Solomon 
[32j[33] and Aliste-Prieto, Coronel, Gambaudo [21 [3], who obtained estimates for 
the rate of convergence to frequency of the number of prototiles per volume for a 
class of domains. They were motivated by questions on bi-Lipschitz equivalence 
and bounded displacement of separated nets, arising from self-similar tilings, to 
the lattice. The reader is referred to remarks at the end of Section 4 for a more 
detailed discussion of these results and how they compare to ours. 

Our goal is a finer analysis of the deviation from the ergodic average, which we 

d-1 

can perform in the case |#2| > 6\ d ■ The main tool here is a family of finitely- 
additive measures associated with the system. It is known that the right and left 
eigenvectors of S corresponding to the dominant eigenvalue 6\ give rise to the 
unique invariant probability measure for the tiling dynamical system. The tiling 
space is locally a product of the "Euclidean leaf" — an open set in M. d — and 
the transversal, which is a Cantor set with the structure of a topological Markov 
chain. The invariant measure is locally the product of the Lebesgue measure on M. d 
and a Markov measure. It turns out that for each eigenvalue 9 of S, that is larger 

d-1 

than 6 1 d in absolute value, one can associate two finitely-additive complex (or 
real signed) measures: one defined on an algebra of sets in M. d including Lipschitz 
domains, and another one defined on the transversal. The latter one yields an 
invariant finitely-additive measure for the dynamical system, if we take a product 
(locally) with the Lebesgue measure. 
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Tilings can be viewed as multi-dimensional analogues of substitution dynamical 
systems. By the Vershik-Livshits theorem [37}I38|. primitive substitution dynam- 
ical systems can be equivalently realized as Vershik's automorphisms correspond- 
ing to Bratteli diagrams. Upper bounds for the deviation of ergodic averages for 
substitution dynamical systems have been obtained by Adamczewski [I]; in the 
related context of interval exchange transformations and translation flows on flat 
surfaces, such upper bounds are due to Zorich [39] and Forni [16]. An asymp- 
totic formula for ergodic integrals for translation flows has been obtained in [11 , 
relying on the construction of a special family of finitely-additive invariant mea- 
sures. In particular, G. Forni's invariant distributions are expressed through the 
finitely-additive measures. Limit theorems for translation flows follow as a corol- 
lary of the asymptotic formula. We should mention that related objects (minimal 
cocycles with a scaling property) for 1-dimensional symbolic substitutions have 
been studied by Kamae and collaborators [HI [21] . 

As we said above, the main difficulty of the multi-dimensional case is due to 
the more complicated behavior at the boundary. In the one-dimensional case, 
finitely-additive measures are directly constructed on "Markovian" arcs and then 
extended to general arcs by exhaustion. In the multi-dimensional case, finitely- 
additive measures are first constructed on tiles, and then the question arises of 
their extension to rectangles, discs and so forth. Note, however, that while the 
boundary of an interval consists of two points, the boundary of a rectangle consists 
of several arcs, and their contribution need not be negligible! 

Our first main result (see Theorem I4.3P is an asymptotic formula for the devi- 
ation of the ergodic average in terms of the finitely-additive measures up to an 
error term, generically of order R^ 1 . Under the additional assumptions that the 
tiles are polyhedral, the similarity map is a pure dilation, and the second eigen- 
value 02 is real, simple, and satisfies 62 > |#3|, we prove that the deviations of 
ergodic averages obey a limit law: more precisely, averages on cubes of side r\ n , 
appropriately normalized, converge in distribution to a non-degenerate random 
variable (see Theorem 16. ip . 

2. Preliminaries 

We begin with tiling preliminaries, following [36], see also |24[ [2B| 130] . We 
emphasize that our tilings are translationally finite, thus excluding the pinwheel 
tiling [27| and its relatives. 
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2.1. Tilings. Fix a set of types (or colors) labeled by {1, . . . ,m}. A tile in R d is 
denned as a pair T = (A, i) where A = supp(T) (the support of T) is a compact 
set in R d which is the closure of its interior, and i = £(T) £ {1, . . . , m} is the 
type of T. (The tiles are not assumed to be homeomorphic to the ball or even 
connected. They may have fractal boundary.) A tiling of W 1 is a set T of tiles 
such that R d = {J{snpp(T) : T G T} and distinct tiles (or rather, their supports) 
have disjoint interiors. 

A patch P is a finite set of tiles with disjoint interiors. The support of a patch 
P is defined by supp(P) = U{ SU PP(^) : T € P}. The diameter of a patch P is 
diam(P) = diam(supp(P)). The translate of a tile T = (A,i) by a vector y 6 W 1 
is T + y = (A + y,i). The translate of a patch P is P + y = {T + y : T £ P}. We 
say that two patches Pi,P2 are translationally equivalent if P<i = Pi + y for some 
g <EM. d . Finite subsets of T are called T-patches. For a set C M d we denote by 

T\n = U{TG T : supp(T) C ft} 

the patch of T-tiles whose supports are contained in Q. 

Definition 2.1. A tiling T has (translational) finite local complexity (FLC) if 
for any R > there are finitely many T-patches of diameter less than R up to 
translation equivalence. 

Definition 2.2. A tiling T is called repetitive if for any patch P C T there is 
R > such that for any x G M. d there is a T -patch P' such that supp(i- >/ ) C Br(x) 
and P' is a translate of P. 

2.2. Tile-substitutions, self-afRne tilings. We study perfect (geometric) sub- 
stitutions, in which a tile is "blown up" by an expanding linear map and then 
subdivided. A linear map : M d — ^ IR^ is expansive if all its eigenvalues lie outside 
the unit circle. 

Definition 2.3. Let A = {Ti, . . . ,T m } be a finite set of tiles in R d such that 
Ti = (Ai,i); we will call them prototiles. Denote by Vj{ the set of patches made 
of tiles each of which is a translate of one of Tj 's. A map to : A — > Vj^ is called a 
tile-substitution with expansion (p if 

(1) supp^Tj)) = 4>Aj forj<m. 

In plain language, every expanded prototile (pTj can be decomposed into a union 
of tiles (which are all translates of the prototiles) with disjoint interiors. 
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The substitution ui is extended to all translates of prototiles by ui{y + Tj) = 
(j)y + uj(Tj), and to patches by w(P) = U{lo(T) : T £ P}. This is well-defined 
due to ([T]). The substitution to also acts on the space of tilings whose tiles are 
translates of those in A. 

To the substitution ui we associate its m x m substitution matrix S, with Sij 
being the number of tiles of type i in the patch ui{Tj). The substitution u is called 
primitive if the substitution matrix is primitive, that is, if there exists k £N such 
that S k has only positive entries. 

Definition 2.4. Given a primitive tile-substitution u), let be the set of all 

tilings whose every patch is a translate of a subpatch of uj n (Tj) for some j < m 
and n £ N. (Of course, one can use a specific j by primitivity.) The set X w is 
called the tiling space corresponding to the substitution. 

Definition 2.5. A repetitive tiling T, such that uj(T) = T for a primitive tile- 
substitution uj, is called a self-affine tiling. The self-affine tiling is self-similar if 
the expansion map of u is a similitude, that is, for some A > 1 we have 

|0O)| = X\x\, for allxe R d . 

The number A is called the real expansion constant, or linear dilatation, of the 
map (p. 

We say that a tile-substitution uj has FLC if for any R > there are finitely 
many subpatches of uj n (Tj) for all j < m, n € N, of diameter less than R, up 
to translation. This obviously implies that all tilings in have FLC, and is 
equivalent to it if the tile-substitution is primitive. 

Remark. A primitive substitution tiling space is not necessarily of finite local 
complexity, see |22[ p. 244] and [T3]. Thus we have to assume FLC explicitly. 
Recently, tiling systems without FLC were studied in [T71 [IS] . 

Lemma 2.6. |26} Prop. 1.2] Let oj be a primitive tile-substitution of finite local 
complexity. Then every tiling S G X w is repetitive. 

2.3. Tile boundaries. For a tiling T denote by 

dT = (J <9(supp(T)) 
TeT 

the union of the boundaries of all tile supports. 
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By the definition of a tiling, dT is nowhere dense in M. d . For self-affine tilings, 
the boundary has zero Lebesgue measure. 

Lemma 2.7. [26], Prop. 1.1] Let T be a self-affine tiling ofR d . Then C d (dT) = 0. 

There is also a kind of "geometric rigidity" : if T is self-affine and dT is piece- 
wise smooth (even piecewise Lipschitz), then it has to be polyhedral. This follows 
from the fact that dT is invariant under the expanding linear map. 

2.4. Tiling topology and tiling dynamical system. We use a tiling metric 
on Xu, which is based on a simple idea: two tilings are close if after a small 
translation they agree on a large ball around the origin. There is more than one 
way to make this precise, and our formal definition is as follows: For 71,75 € X u 
let 

d(T u T 2 ) :=mi{r e (O^" 1 / 2 ) : 3y, \y\<r, supp((7] - y) D 75) D B 1/r (0)}. 
Then c£(7i,75) = min{2 -1 / 2 , d(T, 75)} is a metric on X u . 

Theorem 2.8. [29J (see also [28]). (X u ,d) is a complete metric space. It is 
compact, whenever the space has finite local complexity. The action of M. d by 
translations on X w , given by S i— )■ S — y, y G is continuous. 

This continuous translation action (X w ,R d ) is called the (topological) tiling 
dynamical system associated with the tile-substitution. 

Theorem 2.9. Ifoj is a primitive tiling substitution with FLC, then the dynamical 
system (X^, R rf ) is minimal, that is, for every S 6 X w , the orbit {S — y : y 6 
is dense in X u . 

This follows from Lemma 12.61 and Gottschalk's Theorem [20], see [281 Sec. 5] 
for details. 

Recall that a topological dynamical system is said to be uniquely ergodic if it 
has a unique invariant Borel probability measure. 

Theorem 2.10. If uj is a primitive tiling substitution with FLC, then the dynam- 
ical system (XujMr) is uniquely ergodic. 

This result has appeared in the literature in several slightly different versions. 
We refer to [231 Theorem 4.1] and [28] for the proof. 

Let (j, be the unique invariant measure from Theorem 12.101 The measure- 
preserving tiling dynamical system is denoted by (X w ,M rf ,/i). 
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Lemma 2.11. (see |28} Th. 5.10]) If uj is a primitive tiling substitution with FLC, 

then there exists k G N and T G X^ such that uj k {T) = T ■ 

Combining this with Lemma 12.71 we obtain that C d (dS) = for all S G X w . 

2.5. Substitution action. The substitution uj acts on the entire space X u , and 
it is easy to see from the definition of X^ that uj : X w — > X w is surjective. We 
will address the question of its invertibility, but first record the obvious relation: 



Lemma 2.12. If oj is a primitive substitution with FLC , and jj, is the unique in- 
variant probability measure for the translation action (X^, M. d ), then \x is invariant 
under the substitution action uj (in general, non-invertible) . 

Proof. Consider the probability measure uj*fi = /i o lj~ 1 on X u . It is immediate 
from ([2]) that this measure is invariant under the translation action, hence fi = ui*/j, 
by unique ergodicity. □ 

Definition 2.13. A primitive tile-substitution uj is called non-periodic if all T G 

X^, are non-periodic, that is, T — y = T implies y = 0. If at least one T G X^ is 
non-periodic, then uj is non-periodic by minimality. 

Theorem 2.14 (|35j). The map uj : X w — > X w is injective if and only if uj is a 
non-periodic substitution. 

We assume that uj is non-periodic for the rest of the paper. 

For non-periodic substitutions we have the Z-action generated by uj, along with 
the translation Reaction. It is useful to note that this Z-action is, in some sense, 
hyperbolic, with the orbit of T under the Reaction playing the role of the unstable 
set of T (clear from ([2]) ) , and the the transversal containing T playing the role of 
the stable set, see [1] for details. (The transversal is defined here as the set of all 
tilings which agree with T exactly on some patch, possibly one tile, containing 
the origin in its interior; see Section 4 for precise definitions.) 

The substitution uj can be extended to "super-tiles" and "sub-tiles" of all or- 
ders. More precisely, for any k G Z let 



(2) 



uj(T-y) 



w(T) - 4>y. 



cf> k A := W k m)}T =1 , (<f> k uj)(<p k T t ) = </> h (u{Ti)). 
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This defines a bi-infinite sequence of tile-substitutions and the corresponding 
tiling spaces {X^k^h&L- The subdivision map 

(3) T fc : X^ -)■ X^k-iu 

acts by subdividing each tile according to the rule implicit in the substitution. 
The inverse of the subdivision map is the composition map, which is well-defined 
if and only if the substitution is non-periodic, by Theorem 12.141 For T £ X u 
denote := T, and 

\ T l+ i---T (T), if fc < 0; *" 

Note that is always defined for k < 0, and if oj is non-periodic, then for k > 
as well. 

2.6. Some geometric measure theory. Denote by T~L a the a-dimensional Haus- 
dorff measure (see e.g. [25J for definitions and basic properties) and by C d 
the Lebesgue measure in M. d . A set H C M. d is said to be m-rectifiable for 
m e N, m < d, if U m (H) > and there exist Lipschitz maps hi : R m — > R d , 
i E N, such that 

(4) U m ^H\ Q^(R m )j = 0. 
See e.g. [251 P-204]. 

We will say that an open bounded set O C R d is a Lipschitz domain if there 
exist finitely many Lipschitz maps hi : R d - X -> M d , i < iV, such that 

(5) h*- 1 (an \ (J ^(lOj = 0. 

Thus, the boundary of a Lipschitz domain is (<i — l)-rectifiable. For A C M rf 
denote 

U(A,r) = {x£R d : dist(ar, A) < r}. 
The a-dimensional upper Minkowski content of A is defined by 

M* a (A) = limsup(2r) d - a C d (U( y A,r)). 

It is known [15j that M* m {A) = U m (A) for an m-rectifiable set A. It follows 
that for any (d — l)-rectifiable set A and b > there exists C(A, b) such that 

(6) C d (U(A,r)) < C(A,b)r, for all r G (0, b] . 
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Indeed, we have M* (d - l \A) > 0, hence © holds for some b = bo, and then we 
can simply take 

2.7. Linear algebra. We will need the following well-known result. 

Lemma 2.15. Let J be a Jordan cell of size s with diagonal entries 9, \9\ > 1. 
Then 

' sk s ~ l \Q\ k for k>0; 
s\k\ s - l \6\ k+s fork<0, 
where \\ ■ \\ is the operator matrix norm induced by the Euclidean norm. 

3. Finitely- additive measures on Lipschitz domains. 

Recall that we have the prototiles Tj, i = 1, . . . , m, with supp(Tj) = A4. By 
the definition of self-similar tiling and the substitution matrix, we have 

m 

(7) Aj = + j = l,...,m, 

i=i 

where the union is almost disjoint (up to the boundaries), and Dij are finite sets, 
with Sij = i^Dij. Then we obtain, using Lemma 12.71 that 

m 

£ S ^ £ "(^) = ^(buppCw^)) = C d {\A 3 ) = \ d £ d (Aj), j = 1, . . . , m, 
i=i 

where A is the linear dilatation of (/). It follows that (C d (Aj)) , IL 1 is the Perron- 
Frobenius eigenvector for S*, with the eigenvalue X d , hence the dominant eigen- 
value of S* is 

(8) Q x = \ d . 

Denote by E + the linear span of the Jordan cells of the transpose of the substi- 
tution matrix S* corresponding to eigenvalues greater than 1. Our goal is to define 
finitely-additive (complex) measures for v E E + and T G X u , analogous to 
those from [21 Ell [H], which corresponds to d = 1. 

Recall that for each T £ X w we have a sequence {T (fc) } fce z of "sub-tilings" (for 
k < 0) and "super-tilings" (for k > 0), together with = T, such that 7"( fc_1 ) 
is obtained from by the process of subdivision. They are uniquely defined 
by the assumption of non-periodicity, see Subsection 2.5. 
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Initially, 7- is defined on the following ring of subsets of M. d : 

C+:={J{F = (supp(P)\N 1 )UN 2 : P C T {k) , JVi , N 2 C dT {k) } . 

(Recall that a ring of sets is closed under finite unions and intersections, but not 
necessarily under complements.) In other words, this is the ring generated by 
sub-tiles and super-tiles of T of all orders (referred to as "tiles of order k" for 
k G Z) and arbitrary subsets of their boundaries. The finitely-additive measure 
1S defined on tiles of order k G Z by 

(9) d>+ r (supp(T)) = ((S')H', if 3 A; G Z, y £M d : T = fe (T,) -ye T (fc) , 

and 

$+ r (iV) = ifiVc<9T (fe) . 

Then this set function is extended to finite disjoint unions by additivity. We 
need to show that this definition is consistent; then we get a finitely-additive set 
function on the ring C^-. Clearly, it suffices to verify finite additivity for a tile 
and its decomposition into sub-tiles. We have 

*<r(4? ~ y) = v v where Tj - y G T, 
and hence (J7|) implies 

m m 

E E Kt^~ 1a * + x - v) = Ecsvcs*)- 1 *)* = v i> 

i=l x£Dij i=l 

as desired. Note that (S*) _1 t; is well-defined because v is in the expanding sub- 
space for S*. We showed finite additivity when subdividing a 7~-tile into T^~ 1 ^- 
tiles; the general case follows. 

Observe that v \— > is linear, so we can restrict ourselves to v from a basis 
of E + ; specifically, to a basis of eigenvectors and root vectors of S , associated to 
the canonical Jordan form of S*. 

Definition 3.1. Define the rapidly expanding subspace of E ++ of S* to be 

the linear span of Jordan cells of eigenvalues satisfying the inequality 

\e\ > = x d -\ 

where A is the linear dilatation of (ft, see (Ej). 
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The space E ++ yields finitely-additive measures for which the main contribu- 
tion is in the interior of the set, rather than at the boundary. Heuristically, the 
main contribution of the eigenvalues with absolute values in (1, X d ~ l ) will, on the 
contrary, be concentrated at the boundary. At any rate, it will turn out that the 
contribution of the latter eigenvalues to the deviation of the ergodic average on, 
say, a ball of radius R is bounded above by Ci? d_1 , see ([52]) below. This effect 
was not present in the one-dimensional case, since for d = 1 we have E ++ = E + . 

Let 7 > be the average of A d_1 and the smallest absolute value of an eigen- 
value of S*|£++. Then 7 > A d_1 , and we have 

(10) \\{^) k v\\ < const -7 fe |M|, £;<0, Vu G E ++ , 

where the constant depends only on the matrix S. 

Denote by Q the ring of sets generated by Lipschitz domains in M. d and subsets 
of their boundaries for a fixed d (recall that Lipschitz domains are assumed to 
be bounded by definition). We will show that for v G E ++ the finitely-additive 
measure $^7- can be defined in a natural way on Q. If the tile supports belong to 
Q (in which case they have to be polyhedral, see Subsection 2.3), then this is an 
extension of the corresponding finitely-additive measure defined on C^- But we 
also allow "fractal" boundaries, in which case it is not clear whether the finitely- 
additive measure can be extended to the ring generated by C^- U Q. 

Lemma 3.2. For any v G E ++ , there exist finitely-additive measures defined 
on the ring Q. Moreover, they satisfy the following "cocycle" conditions for any 
n G Q: 

(11) *Jr-y( n ) = Kr( n + v) f° r aU TgXu, ye R d , 

(12) $+ sr (O) = $J w(r) (0(fi)). 

In particular, for an eigenvector v, with S t v = 9v, we get 

(13) KmT)^)) = ^t,T^) = 0*t,T<fl). 

Before the proof, we introduce a construction which will be used throughout 
the paper. It is an efficient hierarchical "packing" of a Lipschitz domain by tiles 
of varying orders. Analogous constructions have been used in [321 133"1 12~1 [51 [3~T]. 

Fix a tiling T G and a Lipschitz domain Recall that T^\u denotes the 
collection of -tiles whose supports lie in f2. Observe that supp(T^|n) G C^- 
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and supp(T (A;+1) |n) C supp(T (fc) |n)- Further, let 

(14) TZ^(Q) := {T E T {k) \n : supp(T) £ supp(T (fc+1) |n)}- 

In words, 7Z^ (O) consists of those tiles of order k in f2 which belong to tiles of 
order k + 1 that do not lie in Q, hence these tiles of order k + 1 must intersect 
the boundary <9f2. Let <i max , <i m i n be the largest and the smallest diameter of a 
T-tile respectively. Since the largest diameter of a 7~( fc+1 )-tile equals d majX \ k+1 , 
it follows that 

supp(ft (fc) (ft)) C U(dn,d max A + ). 

Denote by a mm = mmj< m £ d (Aj) the smallest volume of a prototile. Then the 
volume of an order k tile is at least ammA*, and we obtain 

(15) m (k) (A) < £ d (£/(dtt, d max A fc+1 ))a m i n A- dfc . 
Proof of Lemma \3. 2\ We define for a Lipschitz domain S7: 

(16) *+ r (fi):= lim $J r (supp(r«|n)). 

' fe— >— oo ' 

Let us show that the limit exists (note that we cannot use monotonicity, since 
the values of Q^T need not be positive or even real). By (fl5l) and ([6]), 

#ft( fc >(n) < C(dQ, l)d m axAa m J n A-( d - 1 ) fc = const • A"**" 1 )* 

for k E Z such that c? ma xA A:+1 < 1, that is, for 

fe < - log d max - log A - 1. 

Now, by finite additivity of $^7- on C.£, in view of ©, 

|ci>+ r (supp(r( fc )b))-ci>+ r (supp(r( fc+1 )b))l = |$+ r (su PP (^ fc )(o))| 

< Yl l^r( su PP( T ))l 

re7e( fe )(n) 

< #7^)(Q)-||(S t )Mll 

< const • A^^VlMI, 

for fe < min{0, — logd max — log A — 1}, where we used (fTUj) in the last step. By 
assumption, \*y\ > A d_1 , hence the last expression tends to zero exponentially fast 
as k — > —00, and the existence of the limit in (|16p is verified. We then define 



$J r (n un) = $J r (o) for TV c an. 
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Now let us check finite additivity. If f2i and are Lipschitz domains with 
disjoint interiors such that 

o = Hi u n 2 

for another Lipschitz domain £1, then 

since the former consists of those "F^-tiles which intersect dfli n The 
estimate above shows that the $^1"^— measure of the support of the latter patch 
tends to zero as k — > —00. The definition (|16p then shows 

$+ T (0) = $+ T (0 1 ) + $J T (fi 2 ). 

Thus we can define on a finite union of disjoint Lipschitz domains and subsets 
of their boundaries consistently, and finite additivity follows. 

Formulas (jlip and (|12p are easily verified: they hold on the ring Cj- by defi- 
nition, hence they hold for Lipschitz domains by (|16p . and therefore for all the 
elements of the ring Q. □ 

In the next lemma, we estimate the growth of our finitely-additive measures 
on dilations of a given Lipschitz domain. 

We will often use constants which depend only on the tiling substitution uj 
(which includes all the data for the tiling space X^), and on the domain O. This 
will be often written as C = C(ui, Q) for short. 

Lemma 3.3. Suppose that v £ E ++ , with ||t> || = 1, belongs to the S 1 -invariant 
sub space corresponding to a Jordan block of size s > 1, with an eigenvalue 9 (v is 
an eigenvector if s = 1). Then for a Lipschitz domain £1 and CLr = we have 
for R> 2: 

s-loa „„l 



(17) \^ T (n R )\ ^dilog R) s - L R a , where a x . 

' 1 log A 

for a constant C\ = Ci(ta,tl). 

Proof. Fix a tiling T £ X u . For k E Z we consider patches lZ( k \£l) introduced 
in (|14p . If we fix and let = R£l, then we obtain by (|15p . 



#n( k \n R ) < c d (u(dn R ,d max \ k+1 ))a a l\- dk 

(18) = R*£ d (U(dn,d nm \™R- 1 ))c£*\- dk . 
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Denote 

k R = max{/c e Z : TZ (k) (n R ) / 0}. 
It is clear that d m i n X kR < diam(Q^) = R diam(fi), so 

(19) X kR < Rdiam(n)/d min . 

Thus, in (fTHj) we have d max \ k+l R~ 1 < A diam(fi)d max /d m i n := 6, hence 

(20) #^ (fc) (0 R ) < C" • i? d - 1 A-( d - 1 ) fe , with C = C(dn,b)d m&x Xa^ m . 

where we used ©. Note that the constant b, and hence C , depends only on ui 
and fi. We have by (fTB]) : 

(21) <r(^)= E E <*\M su pp( r ))- 

Therefore, 

(22) |*J r (fi fl )| ^ C"-^" 1 E 



/B ||(s')Mll 

\(d-l)k 
k=—oo 



\k=— oo fc=l / 

(24) < c""-ir*- l |A *' - 1 1 . 

We used ((20]) and © in ([22]), Lemma EH] in ([23]), and the assumption \9\ > A 
in (|24h . Note that the constants C", C" depend only on uj and O (they depend on 
the substitution matrix, which is encoded in uj). We also assumed that k R > 0, 
but this does not lead to loss of generality since it is enough to establish (JTTJ) for 
R sufficiently large. In view of (|19p . we have 



id-l 



\6\ kR log(|9|/A d - 1 ) Log|9j (d 1 v 

— < const ■ R lo § x = const • R lo % x , 



hence the inequality (|24p implies 

|*ir( n «)l ^ ^(logi?)^ 1 ^!^/ 10 ^, 

with Ci = Ci(u;,f2), as desired. □ 

We record the following fact, which easily follows from the proof of Lemma [3.3l 
for future use. The notation T\n R , used in the lemma below, means the collection 
of all T tiles contained in £l R . 
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Lemma 3.4. For a Lipschitz domain Q there exists a constant C 2 = C 2 (u}, Q) > 
such that for all v £ E ++ , with \\v\\ = 1, and for all T £ X u , we have for 

n R = m : 

(25) |^ r (^R) - $+ r (supp(TK)| < C 2 R d ~\ for all R>1. 

Proof. We have 

$+ r (su P p(TK)) =J2 E <r( su PP( T ))- 
Comparing with ()21[) and using (j20|) we obtain, similarly to (|24P : 



l$+ (nB)-*; r (su PP (ria fl )i < 



£ 53 < r (supp(T)) 

fc=-ooTG7e( fc )(n J {) 



- ^(d-l)fc 
fc=— oo 

— 1 k^J 

< C'R^-C" ^ 



fc=— oo 



where 7 > A is from (jTOj) . □ 

3.1. Holder estimates. Next we establish a Holder estimate for our finitely- 
additive measures. Although more general Lipschitz domains could be handled, 
we restrict ourselves to cubes, for simplicity and because the limit law in Section [6] 
below is obtained in this setting. We do not need this result until Section [6l 
Denote 

Q r := [-r/2,r/2] d and A rur2 := Q r2 \ int(Q ri ) for < T\ < r 2 . 
Thus, A ri>r2 is the closed "annulus" between two concentric cubes. 

Lemma 3.5. Suppose that v £ E ++ , with \\v\\ = 1, satisfies S t v = 6v (so that 
9 > A d_1 by the definition of the rapidly expanding subspace E ++ ). Then there 
exists a constant C3 = Cs(uj, £2) > such that for any T G X w and any < n < 
r 2 we have 

(26) |$J r (Q r2 ) - *J r (Qrx)| < C 3 r d 2 ~\r 2 - n)^ 1 ), where a = ^ . 
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Remark. Taking r\ = we obtain the upper bound C^r^, which agrees with 
(fTTl) . since s = 1 for the eigenvector v. 

Proof. We have by finite additivity: 

KriQn) - KriQn) = < r (^i,, 2 ). 
Consider 7Z^ ( A ri jr2 ) as defined in ([H]h recall that 

#K ik) (Ai,r 2 ) < C d (U(dA ri ^,d max X k+1 ))a^ Q \- dk 
by (fl5|) • Clearly, cL4 ri )r2 = <9Qri U dQ r2 . The following claim is elementary. 
Claim. For any r > and i G (0, r), 

(27) £ d (£/(dQ r ,i/2)) < d2 d tr d -\ 

Indeed, we have C d (U(dQ r ,t/2)) < (r + t) d - (r - i) d whence ([57]) follows by a 
simple calculus exercise. 

Therefore, for all < r\ < r2, 

C d (U(dA rur2 ,t/2)) < C d (U(dQ ri ,t/2)) + C d (U(dQ r2 ,t/2)) < d2 d+1 tr^ 1 . 
Thus, 

(28) #K {k \A ri>r2 ) < C{d,uj)\-^ k r d -\ where C(d, u) = d2 d+2 d miiX Xa^ m . 
Let k = max{k G Z : 7^ fe )(,4 ri ,r 2 ) ^ 0}. We have 

*+ r (A^,r a )= E E < r ( su PP( T ))- 

fc=-oo Te72.C*)(^ 1)r2 ) 



Then we obtain, using (|28j) . 

i*; r (A-i,n.)i < c(d, W )r 2 d - 1 e 



\(d-i)k 

fc=— oo 



(29) 



(d-: 

cfauyt 1 ( \Q\ xfc0 



E A(rf -Dfc 

fc=—oo 



1 _ Afzi VA d_1 



It remains to estimate fco- By definition, a tile of order ko must be contained in 
A rij7 . 2 . Let 77 > be such that every T prototile contains a ball of radius r/ in its 
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interior. Then a ball of radius X k °r] must be contained in A ri ^ r2 . It is easy to see 
(we do not attempt to get a sharp estimate here) that 

\ k °T] < T2 — 7*1, 

since the center of the ball must have at least one coordinate in the interval 
[r\ , 7*2] . It follows that 

\0\ \ I T2 — T\ \ ; A 



whence (|29|) yields 

(30) l*J r (^,„)l < ^rf 1 ^ - n )—- {d - l \ 

with the constant C3 depending only on the tiling substitution 00, as desired. □ 

4. Finitely- additive measures on transversals and statement of the 

main theorem 

Recall that the "Euclidean leaf," or the translation orbit, of a tiling T G X w is 
the unstable set for the substitution map oj. The stable leaf is a transversal, which 
we now define, and which is topologically a Cantor set for aperiodic tilings. We 
then proceed to the construction of finitely-additive measures on the transversals. 
This construction is naturally dual to the one in the previous section. 

Definition 4.1. For an admissible patch P of tiles in the space X u the set 

V:={TeI w : PCT} 

is called the transversal associated with the patch P. 

The tiling space X w has a local product structure: 

m 

i=i 

where ~ is a natural homeomorphism, Tj are the prototiles, Aj = supp(Tj), and 
the quotient ~ corresponds to a certain "gluing" along the boundaries of tiles, 
see [1] for details. In fact, X^ can be considered as a translation surface or 
M d -solenoid HUES]. 

If a patch P C T is such that supp(P) contains the origin in its interior, then 
r^p is a stable set of T for oj, in the sense that 

d(u k (T'),UJ k (T)) < c\- k for all T € r w , P) k G N, 



18 



ALEXANDER I. BUFETOV AND BORIS SOLOMYAK 



by the definition of the tiling metric d. 

Now let us derive some properties of the transversals. It is clear that 

(31) w(r U)P ) c r U}UI ( P ), T U)V+P = y + r W) p. 

We don't have equality in the inclusion above, however, we do have 

(32) TfcCrVw^fcp) = r^/s-i^fe-ip, 

where is the subdivision map from ([3|). Then how can we describe (^(r^p) 
precisely? This is the set of tilings T € X u whose "super- tiling" T (1) contains 
the patch <j>P. Thus, we have 

m 

(33) r WiTi = U U MrVJ-x), 

j=l x£D tj 

where Dij are from (|7|), and this is a disjoint union. 

Before we define the finitely-additive measures on the transversals, it is worth- 
while to recall the formula for the unique invariant measure fi. (We know that 
primitive self-affine tiling dynamical systems with finite local complexity are 
uniquely ergodic by Theorem 12.101 ) For a patch P C T £ X w and U C M. d , 
define the set 

Xp t u := {S G X^ : P - y C S for some y G E/}. 

Let 77 > be such that every prototile contains a ball of diameter 77 in its interior. 
It is clear (see e.g. |344 Lemma 1.6]) that the sets Xp^u, with diam(t/) < 77 and U 
open, generate the topology on the tiling space X u . It is proved in [3H Corollary 
3.5] that the unique invariant measure [i satisfies 
(34) 

H(X Pt u) = freq(P) • C d (U) for P C T G X w and 17 Borel, with diam(*7) < 77, 

where freq(P) is the uniform frequency of the patch P in T ■ (The existence of 
uniform frequencies is shown, e.g., in [ 24^ Lemma A. 6].) In particular, we have 

(35) KX Tj ,u) =ireq(T j )-C d (U) 
for a small enough U. It is well-known that 

(36) := (freq(T,)) i < m 

is the Perron- Frobenius eigenvector of the substitution matrix S, normalized by 
the condition (v^\ u") = 1, where = (C d (Aj)) 1 JL 1 is a Perron-Frobenius 
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eigenvector of S*. Here and for the rest of the paper we are using the bilinear 
pairing in C m : 

in 

(37) (v,u) = ^VjUj. 

j=i 

We should also note that there is a notion of transverse measure on a transversal 
r. It is a Borel measure v on B(T) such that v{A) = v(A — y) for every A G B(T) 
and y G M. d such that A — y C T. There is a 1-to-l correspondence between finite 
positive transverse measures and finite invariant measures for the tiling system, 
see Section 5]. In our case this is manifested by (|5^|) and (f33|) . 

Next we proceed to define finitely-additive measures on the transversals T w ^i-x 
for i < m and i£Kf We can also define them on the transversals T Wi p for more 
general patches P, but that will not be necessary. 

For S we have the direct-sum decomposition 

C m = E + ®E~, 

where E + is spanned by Jordan cells of eigenvalues of S with absolute value 
greater than 1. For u G E + , j < m, y G M. d , and k > let 

(38) $-(oo k (T u>Tj _ y )) = (S- k u) r 
We have 

w*(r w , Ti - w ) C T UtTi - x ^=> Ti-xe uj k {Tj - y). 

We claim that for each u G E + , (1381) defines a finitely-additive measure on the 
algebra of subsets of ^T t -x generated by the sets cj^r^Tj-j/), with j < m, k > 
and yeR d such that Tj - x G w* 1 ^ - y). It is enough to verify finite additivity 
in ([33|) . since the general case reduces to it easily. We have 

m m 

(t w ,tJ = «i = s h( s ~ 1u )j = s v®u (w(r w>Tj )) 

i=i i=i 

= E E ^MrVJ-x), 

j = l x£Dij 

as desired. 

We are not going to discuss the extension of <3?^ to a larger class of sets, but we 
will need finitely-additive measures m$-, defined locally as the product &~ x C d , 
on the tiling space X u . In order to make this precise, we define the class of "test 
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functions" which we will be dealing with, and we will define their integrals with 
respect to m*- . 

Definition 4.2. A function f on X u is called cylindrical if it is integrable with 
respect to the unique invariant measure \x and depends only on the tile containing 
the origin, that is, 

3i<m, x£R d , G supp(Tj)-z, T, - x G TnT' =>- f(T) = f(T). 

A cylindrical function may be identified, with a family of functions \jp%\i<irrii 
where fa : Ai = supp(Tj) ->• R, fa G L l (Ai) = L l (Ai,C d ) as follows: 

f(T) = fa(x) ifOeAi-x, Ti-xeT. 

The functions fa are only defined £ d -&.e., which does not cause a problem since 
we will integrate cylindrical functions with respect to <J>~ x C d . The simplest 
cylindrical function is the characteristic function of a prototile T,, which is defined 
by fa = 1, fa = 0, j ^ i. 

Now we define for any cylindrical /: 



(39) m 



m m ,. 

*-(/) :=Y,®u{T^)C d {fa) = J> / fa(y)dy. 
i=i i=i Ja * 

Remarks. 1. Finitely-additive measures m^-, for u G E + , are invariant under 
the dynamics: this follows from (|39p and the fact that 



*«(i/ + r W|P ) = ^(r WlW+P ) = *-(r w ,p). 

2. Let G C m be the Perron-Frobenius eigenvector of the substitution 
matrix S, normalized by the condition (v^\u^) = 1. As already mentioned, 
u^P is the uniform frequency of tiles of type j in the tilings T G X w . Thus, in 
view of (1351) . m ffi - is exactly the invariant probability measure ix on the tiling 
space. For a cylindrical function / we denote by ||/||i its norm in L 1 (X w ,/x); 
observe that 



i=i jA i 



3. Cylindrical functions are not dense in L l (X w , fj,); however, it follows from 
|344 Lemma 1.6] that the set of functions {/ o uo~ k : / is cylindrical, k G N} 
is dense. Thus it is useful to compute m^-(f o uj~ k ) explicitly. If / is the 
characteristic function of a tile Tj, then / o to~ k is the characteristic function of 
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the super-tile 4> k Tj, that is, / o ui~ k (J~) = 1 if and only if T € w fc (T^^. ) — y for 
some y £ cj) k Aj. Thus, it follows from ([38]) that 

™ $ -(/°^) = E( s ^ / i>i°<t>~ k ty)dy = \ dk Y J {s- k u) i fr( x )dx, 

keeping in mind that | det(</>)| = A rf . In particular, if Su = 9u, then 
(40) "»«:(/ °""*) = e-*A*m # -(/). 



Denote by E ++ the rapidly expanding subspace for the matrix S, which is, 
by definitions, the linear span of Jordan cells for S corresponding to eigenvalues 

d-l 

greater than 9 1 d = A In our first main theorem, which we state below, only 
the finitely-additive measures m^- with u in the rapidly expanding subspace play 
a role, since only their contribution dominates the "boundary effects." 

Choose a basis {v ^}Y=i f° r consisting of eigenvectors and root vectors of 
S*, according to the ordering of the eigenvalues 

6\ = \ d > | #2 1 > ••• > \9 m \ 

(the eigenvalues are counted with algebraic multiplicity). We set 

= (C^Aj))™ u 

as discussed above. Then consider the dual basis {u^'}JL 1: so that (v® , u^') = 
5ij. This agrees with the definition of it' 1 ) in (|36p . The vectors {u^'}JL 1 are the 
eigenvectors and root vectors of S, so that Su^ = 9vP> if and only if S t v^' = 9v^' 
(note that we do not need to put complex conjugation, by our definition of the 
pairing (|37l) ). Let t be the dimension of the rapidly expanding subspace E ++ , 
that is, 

\e t \>0^ and |0/ + i|<0^\ 
Then {v^'}j =1 is a basis for E ++ and {u^}| =1 is a basis for E ++ . Denote 

*Jr : =*^), r and *r:=*uw 

Theorem 4.3. Let (X^, M d ) 6e a non-periodic self-similar tiling dynamical system 
of finite local complexity, let \i be the unique invariant probability measure, and 
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let Q be a bounded Lipschitz domain in R d . Then there exists a constant C 
C(oj, O) > 0, such that for any cylindrical function f and any T £ X^: 



(41) 



f f(T-y)dy - C d (n R ) f fdfi-Y, K,r^ ■ m ^ (/) 



n=2 

< CR^HlogRYWfh, for all R> 2, 



where s is the maximal size of the Jordan block corresponding to eigenvalues 

d-l 

satisfying \6\ = 9 1 d (if there are no such eigenvalues, then s = 0). 



Remarks. 1. The second term in (|4ip can be written in a way consistent with 
the sum that follows: for a cylindrical /, 

f fdfi = m 9 -(f) and C d (n R ) = $+ r (ft R ). 



2. We can formally interpret (I4ip also in the case when |^| < &i d ', then I = 1 
and the sum in the formula (1411) is zero. 



It is not hard to extend (|4ip to functions of the form /ow"', where k £ N and 
/ is cylindrical, which form a dense subset of L 1 (X W ,^). In the next corollary, 
for simplicity, we assume that S has no Jordan blocks in the rapidly expanding 
subspace and either Q is the ball centered at the origin, or is a pure dilation. 



Corollary 4.4. Under the assumptions of Theorem \4-3j suppose, in addition, 
that S has no Jordan blocks in E ++ , and the finitely- additive measures $~ 
correspond to eigenvectors o/S* and S respectively, with eigenvalues 6 n , forn < t. 
Moreover, assume that either Q is the ball centered at the origin, or (ft is a pure 
dilation. Then we have for any cylindrical function f and k £ N: 



(42) 



I fou- k (T-y)dy-C d {Sl R )( fou>~ k dn- 



£< r (fifl)-m*-(/°"- fc ) 



n=2 



< CR d - 1 X k (log(\- k R)) s \\f o u 



for all R > 2X k , where C is the constant from Theorem \4-3[ 
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Proof. We have 



/ focj-\T-y)dy = [ f{uj- k T-cp~ k y)dy 

(43) = \ dk f f{cj- k T-x)dx. 

Jij)- k n R 



Observe that 



by the assumption on Q and (/), so we can apply (|4ip . with T replaced by to k T 
and R replaced by \- k R. We have £ d {Q x - kR ) = \- dk C d (n R ), 

by (USD, and 

m,-(/) = ^A- M m $; (/ ^) 
by (|40p . Since everything is multiplied by X dk from ()43j) , all "extra" factors cancel 
out. In the right-hand side of @T} we will get C(A~ fc J R) d ~ 1 ( 1 °g(^ _fc ^)) s |l/lli ! 
which is also multiplied by X dk , and keeping in mind that \i is w-invariant by 
Lemma EJ21 we obtain (g2]). □ 



Next we deduce upper deviation bounds from Theorem 14.31 

Corollary 4.5. Let (X w ,E d ,/i) be a non-periodic self-similar tiling dynamical 
system. Suppose that the substitution matrix S has eigenvalues 9\,...,9 m (real 
and complex), counted with multiplicities and ordered in such a way that 9\ > 
1 02 1 > • • • > \0m\- Further, let s be the size of the largest Jordan block associated 
with the eigenvalues of absolute value \9%\. 

Given a bounded Lipschitz domain Q, there exists a constant C = C(w,f2) > 
such that for any cylindrical function f , with \\f\\\ = 1, any tiling T 6 X u , and 
R > 2 we have 



f(T-y)dy-£ a (fl R ) fdfi 
where a = dlog 1 6*2 1 / log 9\ 6 (d — 1, d). 



CR d -\ i/|0 2 |<Vj 
< \ CR d - l {\ogR) s , if\9 2 \=9^; 
CR a {logR) s -\ if\0 2 \>9^ L , 



Proof. The first two cases are immediate from (|4ip . since then £ = 1. The third 
case also follows from (I4ip , in view of Lemma 13.31 □ 
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It is possible to show that Corollary 14.51 is sharp, at least in the special case 
when the tiles are polyhedral, in the sense that the powers of R in the right-hand 
side in each case cannot be replaced by a smaller power. 

Remarks. 1. There are a number of results related to Corollary 14.51 in the 
literature. When / is assumed to be the characteristic function of a prototile, 
the corollary reduces to estimates of the rate of convergence to frequency for 
prototiles. In the case d = 1 this is essentially the same as estimating symbolic 
discrepancy for substitutions, which was done by Adamczewski pQ. 

Solomon |321 [33] gives deviation estimates similar to ours for the number of 
tiles in a "super-tile" of high order. The tiles are assumed to be bi-Lipschitz 
equivalent to a ball, but the substitution need not be non-periodic. Under these 
assumptions, the estimates are shown to be sharp. 

Aliste-Prieto, Coronel and Gambaudo [2j [3] obtain analogous deviation esti- 
mates. The paper [2], which deals with the d = 2 case, estimates the deviation of 
average from the frequency for general Jordan domains and for very general sub- 
stitution tilings, including non-FLC tilings, the "pinwheel-like" tilings and tiles 

with fractal boundary. However, the extension to d > 2 in [3] handles only the 

i 

case of "small" #2 under the stronger assumption | Qi \ <9^. 

Interest in such estimates was inspired by questions on bi-Lipschitz equivalence 
and bounded displacement of separated nets (also called Delone sets) arising from 
primitive substitutions, like the Penrose tiling, to the lattice in R d , see |12j . 

2. Sadun [31] obtained deviation estimates for the number of patches per 
volume in balls of large radius using rational Cech cohomology, with an error 
term computable from the patterns that appear on the boundary. 

5. Proof of Theorem 14.31 

We will use the following notation: for a set E C W 1 , a tiling T, and a patch P, 
denote by Np(E,T) the number of translated copies of P in the tiling T whose 
support is contained in E. Since T is now fixed, we will just write Np(E) = 
N P {E,T). 

Writing the cylindrical / as a sum over prototiles i < m, we can assume without 
loss of generality that ipj = for j ^ i, and let ifi := ipi. Denote 

!■= [ f(T-y)dy. 
Jn R 
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It follows from the definition of / that if y belongs to a translate of Xi in T, that 
is, y G Ai — x and Tj — x G T for some cc G M d , then 

f(T-y) =ip(y + x), 

and f(T — y) = otherwise. Thus, 

(44) X= ^ / ^(y + x)dy, 

x: (^-x)nti R #« j { jlt - I ) nn * 

where the sum is over x such that Tj — x G T~. Every translate of Tj which is 
contained in contributes £ d (ip) = f A . ^{y) dy to X, and every translate of T< 
which intersects the boundary of contributes at most H^Hi = 
Notice that the number of the translates intersecting the boundary does not 
exceed £ d (U(dQ R , d mSLX ))a^ in . We can write 

c d (u(dn R , d max )) = R d c d (u(dn, d max /R)) < c(dn, ljiUJ^ 1 , for r > d max , 

by ©, hence 

(45) 1 = N Ti (n R )£. d &) + OiR^Wfh), 

where the implied constant in O(-) depends only on Q and u. Thus it suffices to 
prove the desired estimate for Nt^r)- (Note that \C d (ijj)\ < (u\ so 
we will get the factor of ||/||i in the right-hand side of (f4~Ijh ) 

By the definition of the substitution matrix S, we have 
(46) 

w\Tj) - y G T => N T ^ k (A 3 ) -y) = S k (i,j) = (S*)*0'> i) = <(S*) fc e« , e <J)), 

where eW is the standard i-th basis vector. 

Recall that we have chosen a basis {i/ n )}™ =1 for C m , such that {f( n )}^ =1 is 
a basis for the S*-invariant subspace E ++ , and a dual basis {ir ri )}™ =1 . Then we 
have 

m m 



n=l n=l 

Therefore, 

(47) <(S*) fc e« , e^) = ^((S*) V^),, 

n=l 

Next we essentially repeat the construction of Lemma 13.21 and consider the set 
n (k) = -]Z^(n R ) defined by (H). Further, let us write TZ^ = \J? =1 nf \ where 
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71 j is the set of tiles of order k in TZ^ of type j. Let k R = max{/c : ^ 0}. 



We have, in view of (f4U|) and (|47|) . 

A^r) = ^iV Ti (supp(^ fc ))) 



fc=0 



fc=0 j=l n=l 



n=l n=£+lj k=Qj=l 
=: h+h. 

Recall that 

$+ r (supp(r)) = ((S*)M n >)j for T £ T {k) of type j. 
Using this and finite-additivity of ^^-p we can write 

/^^^^(HuppCTIn,)). 

n=l 

By Lemma 13.41 

(49) <D+ r (supp(Th J) = < r ( n «) + f ° r » < 4 

where the implied constant depends only on and oo. Recall that u, 
$-(r W|Ti ). Thus (09]) yields 



(n) 



(50) N Ti (n R )£ d &) = £ < r (^) • m - (/) + J 2 • C d W + 0(R d 



n=l 

with the implied constant that depends only on and u. 
It remains to estimate li- We have 

m kn 

\h\< h^W^Z^^Us^v^l 

n=£+l k=0 

Below we use the notation < to indicate inequality up to a multiplicative constant 
that depends only on and ui. We have 

#n {k) = #n {k) {n R ) < R <i-i\-(d-i)k 

by (jSOD, and 

(51) || (S* ) V n) || < A; 8 " 1 \ {d - 1)k for n > I + 1, k > 0, 
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by the assumption that v^ n \ with n > £ + 1, is in the invariant subspace of S* 
corresponding to eigenvalues 9, \9\ < A 6 * -1 , and s is the maximal size of the Jordan 
block of an eigenvalue 9, \6\ = A 1 . It follows that 

\h\ < R^Y^k*- 1 < R d ~ 1 k s R < R^HlogRY, 

where the last inequality follows from (|19j) . This, together with ()45|) and ()50|) . 
completes the proof of (|41|) in the case when s > 1. If s = 0, that is, all remaining 
eigenvalues are less than A rf_1 in absolute value, then we can replace the right- 
hand side of ([ST]) by 7 fe for some 7 < A d_1 , and use that J2k=o \~ k ( d ~ 1 ^ k X 1 to 
obtain 

(52) \h\<^-\ 

Now the theorem is proved completely. □ 

6. Limit laws for the deviation of ergodic averages 

In order to obtain the limit law, we need to make the following additional 
assumptions: 

(A) the expansion map of the tiling substitution is a pure dilation: 4>(x) = 
Xx, A > 1; 

(B) all the T -prototiles are polyhedral. 

Denote by T the class of bounded cylindrical functions on X u . For any / € T 
and T £ X w , define a continuous function on [0, 1] by 

(53) e n [/,7l(r)= f f(T-y)dy. 

Recall that Q r = [—r/2,r/2] d . We consider r 1— > & n [f, 7~](r) as a random variable 
on (X w ,//) with the values in C[0, 1], endowed with the norm topology. 

Theorem 6.1. Let be a non-periodic self-similar tiling dynamical 

system satisfying the assumptions (A) and (B). Suppose that the substitution 

d-l 

matrix S has a simple positive real second eigenvalue 62 > A^" 1 = 9 1 d , and all 
other eigenvalues are less than O2 in absolute value. Then there is a continuous 
functional f3 : T — > M and a compactly supported non- degenerate measure v on 
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C[0, 1] such that for any f G T satisfying j x f dfi = and (3(f) ^ 0, the sequence 
of random variables 

6n[/,T| 

Km 

converges in distribution to v as n —> oo. 

Remarks. Nondegeneracy of the measure means that if ip G C[0, 1] is distributed 
according to v, then for any ro G (0, 1] the distribution of the real-valued random 
variable (f(ro) is not concentrated at a single point. 

The measure v and the functional (3(f) naturally come from the 2-nd term in 
the formula (|4ip . with Vt = Qi, since the 1-st term in (j41[) is zero. In other words, 

(54) v = the distribution of r h-> f>^ r (Q r ), r G [0, 1], 
ciS ct random variable on and 

(55) /?(/)= m ft -(/). 
Note that 

iw)i<EK (2) i- £d (^)- ii/ii- 

8=1 

by (f39j) . so /3 is a continuous functional on T C L°°(X a; ). 

The theorem in the case d = 1 was established in [10J, and the general scheme 
of our proof is similar. However, it should be emphasized that there are many 
complications because of the "boundary effects" for d > 2. Note that the as- 
sumptions (A) and (B) hold in the one-dimensional case (with connected tiles) 
automatically. 

Proof. We are going to use a basic result (see [3 Th.7.1] or [8]) which says that, 
given a sequence of probability measures on C[0, 1], if their finite-dimensional dis- 
tributions converge and the sequence is tight, then the measures converge weakly, 
which is equivalent to saying that the random variables converge in distribution. 
Recall that a family of probability measures on a separable metric space is tight 
if for every e > there is a compact set such that its complement has measure 
less than e for every measure in the family. 

In view of ([4"Tj) and (fTT|) . we have for / G J 7 , with f x f dfj, = 0, by the 
assumptions on the substitution matrix: 

log 2 



[ f(T-y)dy-$+ T (Q R )- 
■JQr 



m <S>~ (/) 



< CV,ft)iT-l/||i, with a 



log A ' 
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for some 6 £ (0, a) and all R > 2. Therefore, for / £ T with ^ 0, by (15511. 
©n[/,7l(r) *ir(QrAn) 



(56) 



/3(/)^ ^ 



< C(w, n)\~ dn , for all n £ N and r £ [0, 1]. 



Note the following important equality, which follows from (113p and the fact that 
4>{Qr) = AQr = QrA by the assumption (A): 

^MT)iQrx)=02^l T {Qr)- 

Thus, 

*£r(QrA») + 

Observe that r h-> ^2w- n (Tl^ r ) ^ as ^ e distribution of z/ from (|54p for all 
re, since [i is w _1 -invariant by Lemma 12.121 Thus, it follows from (|56p that 
the /c-dimensional distributions of jg^jjgn (©n [/? 7] (n ) , . . . , <5 n [/> 7~]( r fc)) converge 
weakly to the /j-dimensional distributions of (4>J 7 -(Q n ), . . . , 5>2T(Qr fc ))- Further, 
(I26p in Lemma [3.51 shows that the support of v is compact in C[0, 1] by the Arzela- 
Ascoli Theorem. In order to complete the proof, we need to establish (i) tightness; 
(ii) nondegeneracy of the limit measure v. 

6.1. Tightness. The following lemma will imply that the sequence of distribu- 
tions of r i— > 02 n & n [f,T](r) is tight, again by Arzela-Ascoli. In fact, all the 
distributions are supported on a single compact set. 

Lemma 6.2. There exists C{u) and no £ N such that for all f £ T with J f d[i = 
0, for all T £ X U) all n>riQ, and all ri,T2 £ [0, 1], 

(57) \e n[f ,T](r 2 )- e n [f,T](n)\ < c{Xu)moo . , f2 _ rir -(,-D, 

where a = 1 , og e ? . 

log A 

Proof. Let r\ <r<i- We have 

©n[/,T](r 2 )-6 n [/,T](r 1 )= ( f(T - y) dy =: 1. 

JQ\"r 2 \Q\ n r 1 

By the definition of cylindrical functions, there exist tpi £ L°°(Ai), i < m, such 
that 

f(T) = t/)i(x) iff £ mt(Ai -x), 2} - x £ T 
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(we can, of course, ignore the case when x belongs to the boundary of a tile, since 
the boundary has measure zero). Then we have, similarly to 



m . 
1=^2 / ^i{y + x)dy, 

i=l x:(Ai-x)n(Q x n T2 \Q x n n )^0 J(Ai-x)n(Q X n r2 \Q X n ri ) 

where the inside sum is over x such that Tj — x E T ■ For r£N such that A n (r2 
< 1, we estimate X as follows, keeping in mind that ||/||oo = maxj 1 1 "0i 1 1 oo : 

m^^iQx-rAQx^-UWoo = l(\ n r 2 ) d - (\ n n) d ] 

= \ nd {r d -r d )- 
(58) < d\ nd {r 2 - n) ■ 



Observe that 



A' 



\ v log A ' < — ri j 



by the assumption A n < (r2 — ri) , keeping in mind that a 

'x d \ 



ITI / \ d \ 

y < 4/lloc (^J fa - n) < rf||/lloo(r 2 - n)^ 1 ), 

which yields ()57p for such n. 

For n 6 N such that X n {r 2 — ri) > 1, we proceed similarly to the proof of 
Theorem 14.31 and estimate 

m 



i=l 



'■max \\J | oo > 



< C d (U (dQ\n n U «9Q A "r 2 , dmax)) ■ a m in a r 
where a max is the maximal volume of a T prototile. By ([2j 

C d (U(dQ Xnri UdQx« r2 ,d m&x )) < d2 d+1 d raSiX {{\ n r 1 ) d - 1 + {X n r 2 ) d - 1 ) 



We have 



< d2 d+2 d m ^\< d - l \ 



^-n(a-(d-l)) < ^ _ ri \a-(d-l) 



by the assumption A n (r2 — 7"i) > 1- Therefore, 



q—n 
'2 



< const • 



(ra-n)^- 1 ), 
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with the constant depending only on X u , and it remains to estimate 

m 



1=1 



This is done similarly to (parts of) the proof of Theorem 14, 3} with some elements 
from the proof of Lemma 13.51 We proceed to the formal estimate. 

Consider = TZ^iQx^Qx^), and let k = max{k : 11^ ^ 0}. Further, 
let TZj be the collection of tiles of type j in IZ^. For i < m, using ([46 p and 
flUj) , we have, similarly to ([15)1. 



fe 

NrAQ^rAQ^n) = ^iVT l (supp(^ fc ))) 

k=0 

ko m m 

k = j = l 8=1 

(1 m \ ko m 

E+E ^EE^kW 
s =l s=2/ fc=0j=l 

_. rW , 7-0 
— • J l -r J 2 • 



Now, 



I? =u^£ d ( S nMr\ Qxnr2 \Q xnri )), 



hence 



^/«-£ d (^) = 0, 



in view of J / dpt = J2T=i u i ■ C (V'i) = 0. Next, 



1 

(1) /-d/ 



i^i < Eii- (s) iiE#^ (fc) ii( s 

s=2 k=0 
k=0 



by the assumptions on the matrix S 4 , where the constant C4 > depends only 
on the tiling space. We have 

#K {k) < C{d,X UJ ){\ n r 2 ) d - 1 \-^ k < C(d,X a; )A (d - 1)(n - fc) 
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by ([281), hence 

ko—n 



(59) 



i=l 



X d-i 



with a constant C5 > that depends only on the tiling space. Recall that 
r]\ k ° < \ n (r2 —Ti), where r\ is the radius of a ball contained in every T prototile. 
Thus \ k °~ n < r2 ~ Tl , hence the right-hand side of (|59p is bounded above by 

Now, combining everything together, we obtain the desired estimate. □ 

6.2. Nondegeneracy of the limiting measure. It remains to prove that v is 
non-trivial and non-degenerate for every r £ (0, 1]. Assume, to the contrary, that 
for some r we have $^((5,-) = c for /x-a.e. T £ X u . By Fubini, we can find 
T £ X w such that 

(60) £ Q d , Vn £ Z, *J w _„ (r _ ra) (Qr) = c 

Here we use that \i is invariant under translations and under the action of uj~ 1 . 
By (lllj) . we obtain that ( I ) ^ 7 -(Qr + x) = c for all x £ Z d , and then by finite 
additivity, 

^2,T(Qkr) = k2 ®2,T(Qr) = ^ c for ^ N, 

decomposing the larger cube into the union of disjoint translates of Q r . On the 
other hand, 

*J r (QA«r)=^*i w -„ (r) (Qr) = ^C 
by (JT3J) . Now take k = [\ n \ and observe that 

|$+ r (Q A n r ) - $+ r (Q fcr )l < const • A n ( rf - X ) 

by (|26p . This implies that c = 0; otherwise, we get a contradiction for n suffi- 
ciently large, keeping in mind that A d_1 < 62- 

Now suppose c = 0. Then $J fiQk-^r ~ rx ) = f° r ^ ^ N and x £ Q rf by the 
argument as above. Then we can approximate supports of the tiles of T by the 
unions of such cubes to conclude that they also have zero <I>;j~ 7— measure. But this 
is a contradiction, since <&\j-{Ai — y) = Vi, the i-th. component of the eigenvector 
of S* corresponding to 62, if Tj — y £ T . 

Let us explain this more carefully. It is only here that we are using the assump- 
tion that the prototiles are polyhedral. Fix a tile Tj — y £ T and denote by Cl n 
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the union of "grid cubes" 2~ n (Q r — rx), with x S Z, d , whose closure is contained 
in the interior of Ai — y. Then V n := (Ai — y) \ £l n is a Lipschitz domain and 
$J 7 -(O ri ) = by the argument above. We essentially repeat the arguments from 
Lemma 13.31 and Lemma 13.51 and start by writing 

fc 

(61) ^ r (K)= £ Yl ^t( su pp( t ))' 

where k = max{/c : lZ^ k) {V n ) ^ 0}. Next, 

(62) #K (k) (V n ) < C d (U(dV n , d max X k+1 ))a^ n \- dk . 
By construction, int(Ai — y) C U(ft n , 2~ n r^/d), hence 

(63) X ko i] < 2- n rVd, 

where r\ is the diameter of a ball contained in every T prototile, thus d max \ k+1 < 
b\ ■ 2~ n r for k < ko for some b\ independent of n. An elementary argument (see 
Lemma 2.2]) shows that for any union F of lattice cubes in Z d we have 



(64) C d (U{dF,t)) < 2(1 + 2b 1 ) d - 1 tU d - 1 (dF), te(0,h], 

where H. d ~ l {dF) is just the surface area of the boundary. Indeed, for every face of 
OF (say, with the "vertical" normal), consider the "parallelepiped neighborhood" 
of the face, with the vertical side length equal to 2t and the other (d — 1) sides 
of length 1 + 2&x- Clearly, it contains the Euclidean neighborhood of the face of 
radius t for all t < b±, and the inequality (f64"J) follows. Scaling by 2 _ ™r, we obtain 

£ d (U(dft n ,t)) < 2(1 + 2b 1 ) d -Hn d ~ 1 {dQ n ), t £ (0,b 1 ■ 2~ n r\. 

Therefore, for large n, such that d max X ko+1 < 1, we have, in view of ([6]), 

C d (U(dV n , d max A fe+1 )) < C d (U(dAi, d max A fe+1 )) + C d (U(dn n , d max A fc+1 )) 

< C(dA h l)d max X k+l 

+ 2(1 + 2b 1 ) d - 1 d max \ k+1 n d - 1 (dttn). 

It is clear that H^idQ n ) are uniformly bounded in n, since Ai — y is polyhedral, 
and Q n is its approximation by a union of 2 _n r-grid cubes. It follows that 

C d (U(dV n ,d m&x \ k+1 )) < b 2 X k , VkeZ, k<k , 

hence, by ([62]) . 

#n {k \v n ) < b 3 \- {d - 1)k , Vk G Z, k < k . 
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Finally, by (l6T|) and (pj) . 

|*i r (K)l < ^ ^-(d-i)fc^fc 

fc=— oo 

(a \ ko 
jh) <&5(2" n r^r-( d - 1 ). 

Since the latter tends to zero as n — > oo we obtain that 

*+ r (4 - 5 ) = $J r (K) + *J r ( n n) = *Jr(^) = °> 
which is a contradiction. The theorem is proved completely. □ 
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